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Transport coefficients can be obtained from 2-point correlators using 
the Kubo formulae. It has been shown that the full leading order result 
for electrical conductivity and (QCD) shear viscosity is contained in the 
re-summed 2-point function that is obtained from the 3-loop 3PI effective 
action. The theory produces all leading order contributions without the ne- 
cessity for power counting, and in this sense it provides a natural framework 
for the calculation and suggests that one can calculate the next-to-leading 
contribution to transport coefficients from the 4-loop 4PI effective action. 
The integral equations have been derived for shear viscosity for a scalar 
theory with cubic and quartic interactions, with a non-vanishing field ex- 
pectation value. We review these results, and explain how the calculation 
could be done at higher orders. 

PACS numbers: 11.15.Tk, II.IO.Wx, 52.25.Fi 



1. Leading order transport coefficients from the nPI effective 

action 

It is well known that selective resummations can play an important role 
in quantum field theory. A prominent example is the hard thermal loop 
theory, which includes screening effects that regulate infra-red divergences. 
The convergence of such perturbation theories can typically be improved 
using an n-particle irreducible (nPI) effective action. In addition, nPI ap- 
proximation schemes can be used to study far from equilibrium systems. 

Although the nPI effective action is consistent with the global symme- 
tries of the theory, symmetry identities for n-point functions may not be 
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satisfied during intermediate steps of the calculation [U [2] . For scalar the- 
ories, Goldstone's theorem may not be satisfied, and for gauge theories. 
Ward identities can be violated. To address this problem we look at the 
resummed nPI effective action, which is defined with respect to the self 
consistent solutions of the n-point functions. This type of strategy was 
originally proposed by Baym and Kadanoff The resummed effective 
action respects all symmetry properties of the theory, and the n-point func- 
tions which are obtained by functional differentiation satisfy Goldstone's 
theorem, or the Ward identities. 

Transport coefficients can be calculated from the Kubo formulae. Con- 
ductivity and shear viscosity are obtained from the zero frequency limit of 
the 2-point function. This limit produces pinch singularities, which come 
from pairs of retarded and advanced propagators which carry the same mo- 
menta. When integrating a term of the form / dpo G^'^^{P)G°''^'" (P), the 
integration contour is 'pinched' between poles on each side of the real axis, 
and the integral contains a divergence called a 'pinch singularity.' These di- 
vergences are regulated by using resummed propagators which account for 
the finite width of thermal excitations. This procedure introduces extra fac- 
tors of the coupling in the denominators which change the power counting. 
As a consequence, in even the leading order (LO) calculation, there is an 
infinite set of graphs which contain products of pinching pairs that need to 
be resummed. This is accomplished by solving an integral equation whose 
kernel is the square of the matrix elements that correspond to the 2^2 
scattering and production processes. In gauge theories, in addition to pinch 
singularities, the presence of collinear singularities makes 2 — )• 3 scatterings 
as important as 2 — >• 2 scatterings. These collinear terms are resummed by 
another integral equation. The complete LO result is obtained by solving 
the two coupled integral equations. 

The integral equations that resum pinch and collinear singularities can 
be obtained from the equations of motion (eom's) of the 3-loop 3PI effective 
action [ll[5l[6]. Gauge invariance is automatically satisfied, and all LO terms 
appear without the need for any kind of power counting arguments. 

2. Next to leading order transport coefficients 

Little progress has been made on the calculation of transport coefficients 
beyond LO (see however [3 El [9]). Since power counting is notoriously dif- 
ficult, and becomes increasingly complicated at higher orders, nPI effective 
theories could provide a useful method to organise the calculation of trans- 
port coefficients at next-to- leading order (NLO). 

In spite of the complexity of the rePI effective action, a systematic ex- 
pansion can be done in a self-consistent way [lOj. Using the 4-loop 4PI 
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effective action for a scalar theory witli cubic and quartic interactions, with 
non-vanishing field expectation value, one can derive the integral equations 
that would produce the NLO order contribution to the viscosity We 
remark that the numerical solution of these equations is more difficult than 
anything that has been accomplished so far. In addition, renormalizabil- 
ity has not been demonstrated beyond the level of the 2PI effective action 

We use a compactified notation in which a single numerical subscript 
represents all space-time co-ordinates. For example: the field expectation 
value is written := (p{x), the variational propagator is written D12 := 
D{xi,X2), the bare 4-point vertex is written ^1234 := y^{xi,xs,X2,Xi), 
etc. We also use an Einstein convention in which a repeated index implies 
an integration over space-time variables. Using this notation we write the 
classical action: 

Scim = P?2)"']02 - ^C/l°2°3'Al</'203 - ■^V,%(^i(^2<t>3^A ■ (1) 

The 4PI effective action has the form: 

U, V] = S,i [<P]+'-Ty LnD^i (2) 

+ '-Tr 
2 

The terms ^'^[(p, D,U,V] and ^^^^[D,U,V] contain all contributions to the 
effective action which have two or more loops. The first piece ^^[(p, D, U, V] 
includes all terms that contain bare vertices. The diagrams up to 4-loops 
are shown in Fig. [TJ 

There are four eom's which are obtained by functionally differentiating 
with respect to the four functional arguments of the effective action: 

6T[cP,D,U,V] Y ai^nu v\ CX^ 
— = , Aj e {0, D, (7, V\ . (3) 

The equations obtained by varying with respect to {-D, U, V} can be solved 
simultaneously for the self-consistent solutions which are functions of the 
field expectation value: -Dfc/*], C/[</'], y[(t>\- Substituting these self-consistent 
solutions we obtain the resummed action, which depends only on the ex- 
pectation value of the field: 

m = T[^,mM<f\,vm. (4) 

The re-summed propagator and self energy are defined as: 

^(^s*)"' = Ti^f M = im.m-' - ns*] . (5) 



{DU^))-^ [021-01^]] -i^^[c^,D,U,V]-i^^'[D,U,V]. 
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We can derive an expression for the re-summed 2-point function by taking 
derivatives of the resummed effective action and using the chain rule. The 
result has the same form as the Schwinger-Dyson (sd) equation and depends 
on the 3- and 4-point vertices defined in ([6]). We also define a 5-point vertex 
that we will need below. 

n - ^^^^ ^ o -'^^123^ ((^) 

"123 — 7- — ) ^1234 — ^-7 — , ^12345 — — TT • lOj 

d(p3 d(p4 d(p5 

The vertices in equation ([6]) satisfy integral equations obtained by taking 
functional derivatives with respect to the field expectation value of the ap- 
propriate eom ([3]). These integral equations are shown in Fig. [2j One of 
the kernels is shown in Fig. [3l The integral equations for the vertices in ([6]) 
depend on the vertices U and V, which in turn satisfy their own integral 
equations, that are obtained directly from ([3]) and shown in Fig. SI In all 
figures we combine diagrams that correspond to permutations of external 
legs. The full NLO contribution to the shear viscosity can be obtained from 
the Kubo formula where the 2-point function is given by the sd equation 
with vertices obtained by solving the set of coupled integral equations shown 
in Figs. [5] and HI 



3. A new approach to the calculation of the effective action 

It is interesting to consider how the calculation would proceed at higher 
orders. The nPI effective action is defined as the nth Legendre transform 
of the generating functional obtained by coupling to n source terms. A 
direct calculation beyond the 4- loop 4PI level is extremely difficult [14j , but 
it is possible to calculate the effective action without taking any Legendre 
transforms [15j . 

The key to the method is the introduction of a set of fictitious bare 
vertices: to obtain the n-Loop nPI effective action we include in the La- 
grangian the vertices Vj' for j = 3, 4, 5, 6, • ■ ■ n. The inclusion of the non- 
renormalizable interactions (j > 5) is an organizational trick, and these 
vertices are set to zero at the end of the calculation. From this point on, we 
consider (p = ioi simplicity, and use the notation Vj* and Vj to represent 
bare and proper vertices with j legs. Using these fictitious vertices, we can 
show that the eom's and sd equations are equivalent to the order at which 
the truncated theory respects the symmetries of the original theory^ This 

^ The sd equations and eom's are equivalent only when fictitious vertices are used, and 
this result is important only because it allow us to obtain the effective action without 
taking a series of Legendre transforms. It is not true that the non-perturbative 
solutions of a truncated set of sd equations are the same as the solutions of the eom's 
obtained from the nPI effective action. 
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result allows us to construct the n-Loop nPI effective action directly from 
the sd equations. The method has been used to reproduce, with compara- 
tively little effort, the known results for the n-Loop nPI effective action with 
n = 4 and n = 5. The technique has also been used to calculate the 6-Loop 
6PI effective action which is, realistically speaking, impossibly tedious to 
obtain using Legendre transforms. 

The first step is to compare the perturbative expansions of the sd equa- 
tions and the eom's. The complete set of eom's can be written: 

= Vf + fen;. [Vi\ Vk] + fen, [Vu] . (7) 
The functions fcn;[yj'^, V^] and fcnj[Vfc] are defined as: 



i>3: fc„;W»,V..l=i.«-«^, fc,,[Vil=i.«-^ 

where the hats indicate that the 'basketball' diagrams (for example the 2nd- 
5th diagrams in Fig. [1]) which produce the tree terms in ([7]) are dropped. 
The sign difference for the 2-point function and the missing factor 
occurs because of the fact that it is conventional to write the effective action 
as a function of the propagator instead of the inverse propagator. 

We can generate the perturbative expansion of any functional of proper 
vertices by repeatedly substituting ([7D. We can also repackage a pertur- 
bative set of diagrams as skeleton diagrams by repeatedly using the same 
equation in the form: 

= Vj - fen; [Vi\ Vk\ - fen, [Vk] . (9) 

item 1 If we convert a set of skeleton diagram for the vertex Vj into a series 
of perturbative diagrams using ([7]), the leading loop order of the new 
set of diagrams is greater than or equal to the leading loop order of 
the original set. 

item 2 If we include fictitious vertices for 5 < j < n, we can convert 
skeleton diagrams to perturbative diagrams using ([7D, or perturbative 
diagrams to skeleton diagrams using Q, and the leading loop order 
of the new set of diagrams is equal to the leading loop order of the 
original set. 



We illustrate these statements with an example. We use Lpt to indicate 
the loop order of the perturbative expansion. Consider the skeleton diagram 
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shown in part (a) of Fig. \5\ which is of order L = 2. We can expand this 
diagram as a series of perturbative diagrams using equations of the form ([7]) 
which are shown for this example in part (b) of the figure^ The LO term 
is shown in part (c), and is of order Lpt = 2. Thus we have L = Lpt = 2. 
Now consider the result if we set = 0, which means we remove the first 
diagram on the right side of part (62)- In this case the LO term is shown 
in part (d) and is of order Lpt = 3. Thus we see that if the fictitious vertex 
Vg" is set to zero we have Lpt > L. 

We consider truncating the nPI effective action at m-loop orderly] The 
functional derivative of an m-loop graph with respect to the variational 
vertex Vj opens j — 1 loops. This means that an arbitrary m-loop graph 
in the effective action which contains the vertex Vj produces a term with 
C[m,j] loops in the skeleton expansions of the com for the vertex Vj, where 
we define: 

C[m, j] := m — j + 1 . (10) 

Now we consider the Schwinger-Dyson equations, which form an infinite 
hierarchy of coupled non-linear integral equations. They have the form: 

y/rf = y^o + fcnf [v;°, Vk'"^] . (11) 

Although the structure of the sd equations is very different from the eom's, 
it is possible to show |14j : 

item 3 The perturbative expansions of the n-loop nPI eom's and the sd 
equations truncated by setting V^^^ = V^+fc for A: > 1 both match the 
perturbative expansion obtained from the IPI effective action, and 
therefore each other, to order Lpt = C[n,j]. 

This result is perhaps not surprising, since both formalisms should contain 
the correct perturbative physics up to the order at which the calculation 
is done. Item 3 tells us that both the variational and sd vertex functions 
respect crossing symmetry to order Lpt = C[n,j], which we refer to as the 
"truncation order." 

We can formally write equation (jlip as: 

Vf" = + icn'X, Vj:"] + I,[Vi', V,^'] , (12) 

Wi'^ Vk'] ■■= fcnf Vj!''] - Icn'^iVP, Vj:''] . 

^ The propagators in the skeleton diagrams in Figure \E\ also have to be expanded to 
obtain a perturbative diagram. This will produce extra loops that correspond to self 
energy corrections. In this paper we do not introduce notation to distinguish skeleton 
and perturbative propagators in diagrams. 

^ At m-Loops, the nPI effective action is the same as the (n + 1)PI effective action for 
n> m + 1 [TD]. 
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Comparing d?]) and ()12p and using item 3, it is clear that fcnj[Vfc] and 
/j[Vj'^, y^'*'^] must match each other in the perturbative expansion to order 
Lpt = C[n,j]. Therefore we can rewrite Ij[Vi^ ,Vi^'^] as: 

^.■[^°,n'^']=fcn,[<] + extra, (13) 

where the extra term is of order Lpt = C[n,j] + 1. Using item 2 which is 
only true in the presence of fictitious bare vertices, the extra term can be 
rewritten as a series of skeleton diagrams of order L = C[n,j] + 1. Thus we 
have shown that: 

item 4 In the presence of fictitious bare vertices the sd equations can be 
rearranged to have the same form as the nPI eom's, plus additional 
terms of order L = C[n, j] + 1 in the skeleton expansion, which is 
beyond the truncation order. 

It is straightforward to use the result in item 4 to calculate the n-Loop 
nPI effective action from the sd equations. First use ([9]) to remove bare 
vertices in Ij to order jC[n,j] for n > j > 3. Then join the legs of the 
resulting equation for /j, set = for n > / > j ' + 1, and multiply each 
diagram by the symmetry factor: 

S[j] = s{l/vj){l/jl), (14) 

where s is the numerical factor in front of the diagram in the Ij , and vj — 1 
is the number of times the vertex Vj appears in the diagram. 

It is clear that joining the legs in diagrams in the eom's will produce 
all graphs in the effective action at a given order. The trick is to obtain 
the correct symmetry factor. Consider starting from a known result for the 
effective action <I>'°*, taking derivatives of each graph with respect to each 
variational vertex, and trying to reconstruct the effective action by joining 
the legs in each of these eom's. A given graph in the effective action will 
produce contributions to the eom's of each vertex it contains. In order to 
produce the correct symmetry factor when joining legs, we must drop the 
corresponding contribution in all but one eom, which we can take to be the 
eom for the largest vertex present. This is accomplished by imposing the 
condition Vi = for n > I > j + 1. If the largest vertex appears in a given 
diagram in the effective action more than once, the graph that is produced 
by joining the legs will have a symmetry factor that is too large by a factor 
equal to the number of times the vertex appears. The correct symmetry 
factor is produced by including the factor (l/vj) in equation ([HI). 

Following the simple procedure described above, one can reproduce the 
known results for the 4-Loop 4PI and 5-Loop 5PI effective actions. The 
method has also been used to calculate the 6-Loop 6PI effective action |15| . 
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4. Conclusions 

The LO calculation of transport coefficients involves infinite resumma- 
tions of diagrams that can be identified by power counting arguments. The 
integral equations that resum these diagrams are produced directly from 
the equations of motion of the 3-loop 3PI effective action. This result leads 
to the conclusion that nPI effective theories provide a natural framework 
for the calculation of transport coefficients beyond leading order. The cal- 
culation of next-to-leading order shear viscosity in scalar theory can be for- 
mulated using the 4-loop 4PI effective action. Higher order effective actions 
can be calculated in a efficient way by working with the Schwinger-Dyson 
equations. 
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Fig. 1. 4-loop 4PI effective action. 





Fig. 2. The integral equations for the vertices £7, 5* and Q defined in equation ([5]). 
The kernels are drawn as grey boxes. 
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Fig. 3. One of the kernels in Fig. [51 
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Fig. 4. Integral equations for the variational vertices U and V from the 4-Loop 4PI 
effective action. 




Fig. 5. Diagrams used to explain items 1 and 2 in section [31 



